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Superconducting Gap Modulation in Weak Stripe States 
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The superconducting gap modulation is investigated in the presence of a weak stripe struc- 
ture, using the Bogoliubov-de Gennes theory on the two-dimensional Hubbard model with 
nearest-neighbor site pairing interaction. We calculate the local density of states and discuss 
the recently observed scanning tunneling spectroscopy spectra with four lattice periodicity on 
Bi2Sr2CaCu208+5. We also consider the spectral weight in the reciprocal space, where the 
Fermi surface and the superconducting gap are modulated by the band folding effect of the 
stripe structure. 
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Recently, much attention has been focused on the 
stripe state of underdoped high-Tc cuprates. The stripe 
state was proposed to explain the static magnetic in- 
commensurate structure observed in elastic neutron scat- 
tering experiments on La2-2:Sr2:Cu04 (LSCO)^""^^ and 
Lai.6-£cNdo.4Sra;Cu04.^'^^ It is considered that doped 
holes are localized in the stripe region, which contributes 
to one-dimensional (ID) metallic conduction,^) and the 
outside region of the stripe is an antiferromagnetic (AF) 
insulator. In YBa2Cu307_5, inelastic neutron scatter- 
ing experiments reported incommensurate fluctuations, 
which is consistent with the stripe concept. 

Recently, the incommensurate charge modulation was 
also reported in Bi2Sr2CaCu208+5 (BSCCO). Scanning 
tunneling spectroscopy (STS) revealed modulated quasi- 
particle states with four unit-lattice periodicity sur- 
rounding a vortex core,^*^-* and recent STS experiments 
reported that the four unit-lattice periodic charge mod- 
ulation survives even at zero field without a vortex. -'^^^ 
In the Fourier component of the checker-board pattern, 
the peak intensitiy at either 7r(±-j, 0) or 7r(0, ±i) is dom- 
inant (the lattice constant is taken as unity throughout 
this paper), meaning that we can assume a vertical stripe 
state and that the x- and the y-direction stripe layers are 
stacked layer by layer. 

Our study is based on the self-consistent Hartree-Fock 
(HF) theory of the Hubbard model. It is believed that 
the stripe concept is valid beyond the HF approxima- 
tion.^^) We can consider the metallic stripe state under 
the HF theory if we consider the realistic Fermi surface 
topology. ^''^ -^^^ Once the density of states (DOS) remains 
at the Fermi energy (i.e., metallic state), we can pro- 
duce superconductivity by introducing the pairing inter- 
action, at least as a phenomenological model. In 
the previous study on the superconducting state and the 
vortex state in the stripe state, we considered the case 
where the stripe order is stronger than the superconduc- 
tivity."'^^' However, experimental results on BSCCO 
suggest that the stripe order is not so great. Using the 
same method as ours, the stripe modulation around the 
vortex or the impurity center has been reproduced for 
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the weaker stripe case.^°) Therefore, it is valuable to 
study the STS spectrum for the weak stripe state based 
on this theory. 

In this work, we investigate how the superconducting 
gap is modulated by the stripe order with four (eight) 
unit-lattice charge (spin) periodicity at zero field. We 
use the Fermi surface topology and the electron filling 
~ 0.88 appropriate to optimally doped BSCCO, and 
study the weak stripe state in comparison with the su- 
perconductivity. We calculate the local density of states 
(LDOS), and discuss the STS spectrum on BSCCO. We 
also analyze the spectral weight, which is observed by 
angular- resolved photoemission spectroscopy (ARPES), 
to discuss how the Fermi surface and the superconduct- 
ing gap in the reciprocal space are modulated by the 
weak stripe structure. 

To consider the weak stripe modulation, we begin with 
the conventional Hubbard model on a two-dimensional 
square lattice, and introduce the mean field rii.o- — 
{a\ ^ai_a) at the i-site, where cr is a spin index and 
* — (*x,ij/)- We assume a singlet pairing interaction 
gs between nearest-neighbor (NN) sites, because the co- 
herence length is short in high-Tc cuprates. This type 
of pairing interaction produces d-wave superconductiv- 
ity with the pairing function cos — cos ky . This gap 
function is confirmed by the ARPES experiment. The 
decomposition to the singlet and the triplet components 
of the pairing interaction is explained in ref. 23. Thus, 
the HF Hamiltonian is given by 



+ 9s^{^,^aJ,ia,^] + Aj- ja]_^aj J, (1) 

with a creation (annihilation) operator {ai,a)- For 
the transfer between nearest, second and third neigh- 
bor pairs (j,j), iij = t, t' and t" , respectively. Then, 
the dispersion is given by e(fc) — ~2t{coskx + cos ky) — 
At' cos kx cos ky — 2t"{cos2kx + cos2fcj^). The singlet 
pairing potential is defined as A^-^i = gs{{aj^iai^i^) — 
(aj^ai.i)). When rii^^ and Aej(= Ajj+g) have the A^- 
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site periodicity along the y-direction, we can write 



n 



0<1<N 



'na,lQ, 

^e,j= ^"^'"'^e^iQ {e = x, y), 



0<1<N 



1 

2' 2 



(2) 
(3) 

(1,0) 



with the ordering vector Q = 2'k{ 
and y = (0, 1), respectively mean the next sites in the x- 
and y-directions. There is a relation A_gj = in 
the singlet pairing. 

After the Bogoliubov transformations 




7'^) and at = J2e(^^,rU + u'ljl'l), 



where e is the label of the eigen states, we use the 
Fourier transformations from Uej {ve,j) to u^^k {ve,k)- 
We write fc = fco + inQ {m = 0,1, ■ ■ ■ , N — 1), where 
fco is within the reduced Brillouin zone of size (2^)"^ /N. 
Then, the Bogoliubov-de Gennes (BdG) equation in re- 
ciprocal space is given by 



E 



= E, 



K 



D, 



t 



-K 



Fig. 1. (a) Profile of the charge density, rn at the (0, ij/)-site. 
U = 4.2t (solid line) and 4.3t (dashed line), (b) Spin density 
{—ly^+^v Si. (c) Superconducting pair potential for the d-wave 
^d,i (circle), the s-wave A^^j (square) and the p^-wave Ap^j (tri- 
angle) components. Dotted line shows Aq for the uniform state 
without stripe {U = 0). 



spin electrons, respectively. Then, the LDOS 
N{r,E) = -'K-Hmguir,r,iuJn E + iO+) + 

'jT~Hmg22{r, r, — ia;„ E + iO+) is reduced to 



i,m,m' 



N{r,E) = J2e'''^'^N{lQ,E), 



(4) 



where K^^m,m' = [e{kQ + mQ)- lj]6m,m' +Un_^^(^rn-m')Q 

with chemical potential /i. Dm,m' = 9sJ2e=x,yi^'''''°~^"^'^^'^+ 

(,-i(feo+mQ)-e)^.^^^_^,^^ with Ue,k = Ua,ko,m and Ve,k = 

Va,ko,m- The eigcn-statc e is labeled by feo and the eigen 
state a of eq. (4). Then, the wave function of the eigen 
states is given by 



with the Fourier component 

N{IQ,E) = ^ J2 K,ko,Ako,m+l 



S{E- 



ko,a,m 

+ K,k„-m'>J'^'t'o,m+lS{E + Ea,ko)] 



(10) 

Ea,ko) 
■ (11) 



Similarly, from the Green's function g{k,iiOn) = 
Nk~^ X)r r' e^''^ *^''"'' ^g{r, r', iLOn) in the fc-space, we ob- 



Ua,koir) 



E 



A(ka+mQ)-r 



tain the spectral weight as^^^ 



-1/2 



E^«' 



fco,m^ 



i(fco+mQ)-r 



(5) 
(6) 



A{k,E) 



■ E[ 



^5{E — Ea,ko) 



a^ko^m I 

+ K,ko,mf'S{E + E^^koMko +mQ- k). (12) 



where Nk = J2i 1 = Sfco m 1- obtain the same 



BdG equation as eq. (4) for 

but with the eigen-energy 
/* 

^ a,ko,m^ 



a,ko,'m 



and 



Then, 



/* 

^ a,fco,m 
I '* 

^ a,ko,m^ 



-E'^ feo ^' '^^^y respectively, correspond to 
''a,ko,m, Va,ko,m,, Ea,ko and 7. When there appears spin 

The 



order, u ^ u' , v ^ v' , and E ^ E', generally, 
selfconsistent conditions are as follows, 



Ae.iQ = — 



1 

Nk 



E 

a,ko,m 



[f{E^,ko) e'C^^+^^+'^Q) 
/(-i^^.feJe-'C'^+^Q) ' 



1 
1 



E 

a, fco ,m 



'a,ko,m''^a,ko,m+h 
+lf{Ea,ko), 



(7) 
(8) 



ni,lQ = ]^ E '"a,ko,m'"a,ko,m+lf{-Ea,ko)^ (9) 



with the Fermi distribution function f{E). 

The thermal Green's functions are given by gii{r,r', ioj 

Efco,a "a,fco (^X,feo ir')/iiuJn " £^a,fco ) and 522 [v, v' , l^n) = 

T,ko,a Vc.M{r)vl^ko(r')/{iuJn - E^^ko) for up- and down- 



To reproduce the Fermi surface topology of BSCCO, 

we set t' = -OMt, t" ^ 0.23t and fi OM. The 

chemical potential /k is tuned in order to set the spatially 
averaged electron density ?ii ~ 0.88. The essential results 
of this study do not significantly depend on the choice of 
these parameter values. We consider the case of a vertical 
stripe state with the spin periodicity iV = 8. The stripe 
states appear for U > Uc ^ 4:t on the BSCCO-type Fermi 
surface when g^ = 0. The critical Uc is larger than that 
of the LSCO-type Fermi surface.^'*' ^^-^ We consider the 
superconducting state with = O.St and T ^ 0. In 
this case, the stripe state appears for U > 4. If. We 
report the result for U = A.2t and 4.3i. Since the stripe 
modulation is increased with raising U , the effect of the 
stripe is eminent for U = A.2>t. 

Figures 1(a) and 1(b), respectively, show the profiles 
of the charge density rii = + n,;| and the spin den- 
sity Si = ^[n-i] — nil). The domain wall of the AF is 
located at the bonds between ij,=3 and 4, and between 
iy—7 and 8. Since the bond-centered stripe has lower 
_^ energy than the site-centered stripe in our parameters, 
we report on the former case. Qualitatively, the same 
results are also obtained for the site-centered stripe case. 
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(a) (b) 

Fig. 2. Local density of states N{ri,E) at tlic stripe site iy = A 
(solid line) and the outside site iy = 5 (dashed line). We also 
show the 2Q-Fourier component iV(2Q, £;)e'5"/* (dotted line). 
U = 4.2t (a) and 4.3t (b). 



At the stripe sites next to the domain wall, 1 — is 
larger, i.e., doped holes are weakly accumulated. Cor- 
responding to this stripe structure, the pair potential 
of the superconductivity is also modulated as shown in 
Fig. 1(c), where we show the components of d-wave 
Ad^i = (As,j -I- A_a,i - Ay^i - A_ji,i)/4, extended s- 
wavc As J = {Ax.i + A^x,i + Ay_i + A^y^i)/A andpy-wave 
Ap^i = (Ay^i — A_y_i)/2. The superconductivity is sup- 
pressed in total by the stripe formation. At the hole-rich 
stripe site, A^.^ has large amplitude. When Ad.i has a 
spatial variation, small Ag^i and Ap^j are also induced. 

Figure 2 shows the LDOS at the stripe site and the out- 
side site. In the presence of the stripe structure, there ap- 
pear two peaks inside the gap of the d-wave superconduc- 
tivity. These low energy peaks become larger at the out- 
side site with weak A,;.,;. At the stripe site, the peak at 
the gap edge becomes sharp. Wc also present the Fotiricr 
component N{2Q, E)e'^''/'^ oc N{iy = 5,E) - N{iy = 
4:,E), showing the difference of the LDOS between the 
stripe site and the outside site. The factor e'""^"/^ is 
multiplied so that the Fourier component becomes real. 
Recently, ref. 25 reported the calculation of the 2Q- 
componcnt assuming carge density wave. Since the stripe 
site gains a large normal state DOS with hole accumula- 
tion, the DOS above the superconducting gap is larger at 
the stripe site in Fig. 2. When the calculated spectrum is 
compared with the STS spectrum, they both show 
that the 2Q-Fourier component is positive (negative) at 
lower (higher) energies. However, the weights of the pos- 
itive and negative parts of the 2(5-component are quan- 
titatively different. In the STS spectrum, the weight of 
the negative 2Q-component at higher energies is small, 
because the peak of the gap edge is broad at the site 
with a large superconducting gap. In our calculation, 
the peak at the gap edge becomes sharp, and the weight 
of the negative 2Q-component is large. We give some 
discussion as follows with respect to this inconsistency. 

(i) In order to reproduce the observed STS spectrum, 
we need to consider additional effects of broadening the 
sharp gap edge at the site with a larger superconduct- 
ing gap, for example, broadening by impurity scattering, 
whose contribution may be larger at the site with a large 
normal state DOS. 

(ii) The peak inside the gap (the peak at gap-edge) at 
the outside (stripe) site also appears at the stripe (out- 
side) site, meaning that the peak state extends to the 




(b) kx/n 



Fig. 3. Spectral weight at Fermi energy in the regions | hx \ ^ tt 
and \ky\ < n. We plot integrated A{k,E) within the energy 
range |E - /i| < 0.5t. U = 4.2t (a) and 4.3t (b). (c) Origi- 
nal Fermi surface (bold line) and ±Q-shifted ones (dotted and 
dashed lines). 



adjacent outside (stripe) site. However, we see only one 
peak (higher or lower energy) at each site in the STS 
spectrum. This suggests that the coherence of the peak 
state is restricted to within a few sites by the smearing 
effect of (i), or by other effect which confines the peak 
state to within a few sites. 

(iii) We assume that the pairing interaction is con- 
stant. Then, the site dependence of A^^j comes from 
the site dependence of the DOS. We can consider the 
possibility that the pairing interaction depends on the 
local hole density, having an additional effect on the gap 
inhomogeneity. This may be related to the inter-site co- 
herence of the peak state as noted in (ii) . 

(iv) In our calculation, we assume the incommensurate 
spin order in addition to the charge order. The STS 
observes the charge order with the ordering vector 2Q. 
However, the spin order with the ordering vector Q has 
not been confirmed experimentally. We expect that the 
ordering with Q can be detected by ARPES, as discussed 
later. 

Next, we discuss how the gap modulation due to the 
stripe structure appears in the reciprocal space. The 
wave number fc-resolved DOS is the spectral weight 
A{k, E), which is obtained with the same wave function 
used in the LDOS calculation. When the system has 
a periodic structure, the band folding effect appears in 
A{k, E). Figures 3(a) and 3(b) show the spectral weight 
at the Fermi energy /i in the stripe state, integrating 
A{k, E) in the energy width ±0.5i around fjt. In addition 
to the DOS at the original Fermi surface, there appears 
small DOS near (tt, 0) and (0, tt) in the stripe state. As 
schematically presented in Fig. 3(c), this is due to ±Q- 
shifted Fermi surfaces. To induce new DOS by band 
folding, finite DOS is necessary near the Fermi energy 
in the original dispersion e{k). Around (7r,0) and (0,7r), 
there is large DOS at £^ ~ — l.Ot due to the van Hove 
singularity in our dispersion. It is worth noting that the 
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Fig. 4. The superconducting gap along the Fermi surface. We plot 
A{k, E) at fe along the Fermi surface. Wave number fe is specified 
by the angle 6 from (7r,7r), as shown in Fig. 3(c). U = 4.2t (a) 
and 4.3t (b). 



Fermi surface near (tt, 0) and (0, tt) is still the subject of 

discussion, because additional spectral weight appears 
there. The spectral weight for the (tt, 7r)-shifted 
Fermi surface could appear due to AF fluctuations as 
observed in the ARPES experiment. If the AF fluctu- 
ation is changed to an incommensurate one with ordering 
vector Q, we expect to observe the spectral weight for a 
Q-shifted Fermi surface instead of a (tt, 7r)-shiftcd one. 

To discuss the superconducting gap structure in the 
reciprocal space, we show A{k, E) along the Fermi sur- 
face in Fig. 4. The DOS vanishes when — /i| is smaller 
than the superconducting gap, and has a peak at the gap 
edge. In the uniform state without the stripe, the gap is 
given by | coskx — cosfcj,|, which has a node at 6* = 45°. 
In the presence of the stripe, the gap is modulated as 
shown in Fig. 4. Around 9 = 15°, the Fermi surface DOS 
vanishes upon the reconnection between the original dis- 
persion and the ±Q-shifted dispersions. For 9 > 60°, 
^-dependence of the gap becomes flat. When the x- and 
the ^/-direction stripe layers are stacked layer by layer 
in the sample, or when there is a domain structure for 
the X- and the y-direction stripe regions within a layer, 
we observe an overlap of the spectral weights for both 



direction stripe states, i.e., [A{k, E)+A{kx ^ ky, E)]/2. 

When the stripe modulation becomes strong with in- 
creasing U in our system, the DOS at the Fermi energy 
is reconstructed and they form a ID Fermi surface at 
kx = ±7i'/4, but the DOS vanishes near (7r/2,7r/2), as 
suggested in rcfs. 15 and 29. In this case, the supercon- 
ducting gap in the LDOS N{ri, E) has a small full gap, 
because the low-energy state near the gap node around 
(7r/2,7r/2) vanishes in reciprocal space. This spectrum 
is seen in Fig. 3 of rcf. 17. There, the spectrum at the 
site far from the vortex is almost the same as that of the 
zero-field case without a vortex. 

In summary, we investigate the superconducting gap 
modulation in the presence of a weak stripe structure, 
using the BdG theory. From the same wave functions, 
the LDOS related to the STS spectrum and the spec- 
tral weight related to ARPES are calculated. In the for- 
mer, there appear new peaks inside the superconducting 



gap. In the latter, there appears the eff'ect of band fold- 
ing caused by the stripe structure. On the basis of our 
results, we have discussed the STS spectrum with four 

unit-lattice periodicity on BSCCO.^^'^^^ 
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